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Abstract 

Coplanar-grid (CPG) particle detectors were designed for materials such as CdZnTe (CZT) in which charge 
carriers of only one sign have acceptable transport properties. The presence of two independent anode sig- 
nals allows for a reconstruction of deposited energy based on the difference between the two signals, and a 
reconstruction of the interaction depth based on the ratio of the amplitudes of the sum and difference of the sig- 
nals. Energy resolution is greatly improved by modifying the difference signal with an empirically determined 
weighting factor to correct for the effects of electron trapping. In this paper is introduced a modified interaction 
depth reconstruction formula which corrects for electron trapping utilizing the same weighting factor used for 
energy reconstruction. The improvement of this depth reconstruction over simpler formulas is demonstrated. 
Further corrections due to the contribution of hole transport to the signals are discussed. 
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1. Introduction 

CZT as a 7-ray detector material has a number 
of attractive properties such as a large band-gap 
and good performance at room temperature. The 
COBRA experiment [ 1 ] uses CZT detectors in a 
search for neutrinoless double-beta decay (Oz//3/3) 
due to the presence of several candidate iso- 

topes in CdZnTe and its low natural background ra- 
dioactivity. Currently COBRA utilizes a 4 x 4 x 2 
array of 1 cm 3 CPG detectors operating under low 
background conditions at the Gran Sasso National 
Laboratory (LNGS), with plans to upgrade to 64 de- 
tectors in early 2013. 

Good energy resolution is an important consid- 
eration in gamma spectroscopy. It is also key for 
COBRA in order to enhance the sensitivity to a 
0vf3/3 line (produced by the summed energy of the 



•Corresponding author. Tel 49 351 46334568; fax 49 351 
46337292. 

Email addresses: 
matthew_christopher . f ritts@tu-dresden . de 
(Matthew Fritts), 

juergen . durst@physik . uni-erlangen . de (Jiirgen 
Durst), thomas . goepf ert@tu-dresden . de (Thomas 
Gopfert), Thomas . Wester@physik . tu-dresden . de 
(Thomas Wester), zuber@physik . tu-dresden . de (Kai 
Zuber) 



two betas emitted) and to distinguish it from other 
backgrounds, including the continuous spectrum 
of the competing process of neutrino-accompanied 
double-beta decay (Ivfifi). The depth information 
provided by the CPG design is also very important, 
for two reasons. First, an important category of 
background is a- or /^-radiation at the cathode and 
anode surfaces. The short penetration of such back- 
ground allows for its efficient discrimination, if the 
interaction depth is well known. Second, for certain 
regions in the detector distortion of the energy recon- 
struction is an unavoidable result of the CPG design. 
Depth information can be used to identify events in 
these regions and remove them from consideration 
as Oj//3/3 candidates. For both of these phenomena 
analysis cuts are necessary, and an accurate assess- 
ment of the selection efficiency of such cuts requires 
the most accurate depth information achievable. 

Such considerations give strong motivation for an 
accurate depth reconstruction for CPG detectors for 
their use in the COBRA experiment, and for all ap- 
plications in which accurate knowledge of the ac- 
tive volume is required. In addition, due to the large 
number of detectors used in COBRA, an automatic 
method for generating depth formulas unique to each 
detector is highly desirable. Material properties and 
optimal operational parameters differ significantly 
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from detector to detector, necessitating unique for- 
mulas for both energy and depth reconstruction. 



2. CPG principles 

The CPG detector design was introduced by 
P. Luke in 1994 [2]. A proper understanding of the 
design principles begins with a treatment based on 
the Shockley-Ramo theorem. For a more detailed 
description of the application of the Shockley-Ramo 
theorem to CPG detectors see Q . 

In a CPG detector such as those used in COBRA, 
the cathode is patterned on one side as a uniform 
rectangle, while on the opposite side two isolated 
anode grids are formed in the shape of interlocking 
combs. Figure Q] is a schematic of a COBRA CPG 
detector, manufactured by EI Detection & Imaging 
Systems ||4J]. In operation a large (approximately 
1 kV) bias (referred to as HV) is applied between 
the cathode and one anode, while a smaller (typi- 
cally 50-100 V) bias (GB) is applied between the 
anode grids. The anode held at higher potential is 
referred to as the collecting anode (CA); the other 
anode is known as the non-collecting anode (NCA). 
Electrons excited into the conduction band by a par- 
ticle interaction somewhere in the bulk of the de- 
tector drift straight toward the anode plane (in the 
negative- z direction) until they get close to the an- 
ode grid rails, at which time they are diverted by the 
near-anode field to be collected in the CA. 

A first step in the application of the Shockley- 
Ramo theorem is to calculate the weighting poten- 
tials yl p. 813-8]. Figure [2] shows the weighting po- 
tentials of the CA and NCA along a plane through 
the center of the detector, calculated using the elec- 
trode geometry of COBRA detectors Q At most lo- 
cations in the detector, the two weighting potentials 
are nearly equal, rising with a slope of 1/2 from the 
cathode plane to the anode plane. At the position of 
a CA rail (normally the final location of the mobile 
electrons created by ionization), the CA and NCA 
weighting potentials are 1 and 0, respectively, by 
definition. 

These properties of the weighting potential allow 
for a difference signal based on the raw CA and 
NCA signals that is proportional to the charge from 
ionization and independent of the interaction depth. 
From the Shockley-Ramo theorem one can calculate 




Figure 1: Schematic electrode layout for CPG de- 
tectors used by COBRA. Two interleaved anodes 
and an outer guard ring (uninstmmented) are pat- 
terned on one side (top); the other side (bottom) is 
a uniform cathode. Shown here is the convention 
used in this paper for the z axis (interaction depth). 
The dashed lines indicate the section over which the 
calculated weighting potential is shown in Figure [2] 



1 In current operation the guard ring electrode (black in Fig- 
ure [T} is left unconnected. This information is incorporated into 
the weighting potential calculations. 
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Figure 2: Weighting potentials of CA and NCA as 
calculated for the COBRA CPG detector geometry 
across a plane through the center of the detector. 
The coordinates are normalized to the detector di- 
mensions, thus ranging from to 1 . 



the amplitudes of the two anode signals produced by 
electrons: 

Aqca = ~Qo(*o + i) (l) 

AqNCA = 7jQo(zo - l) (2) 

Here Aq represents the change in induced charge, 
and thus the amplitude of the signal, for the corre- 
sponding anode. Qo is the magnitude of the mobile 
charge produced by the interaction, zq is the distance 
between the interaction location and the anode plane 
expressed as a dimensionless fraction of the detector 
length, zo is referred to as the "interaction depth" (or 
simply "depth"). 

The depth dependence in the two raw anode sig- 
nals is removed by subtraction: 

Aq C A - Aq NC A = Qo (3) 

Note that Aq^cA is always negative, so the equation 
represents a sum of the absolute amplitudes. By cal- 
ibrating this signal using sources of known spectra 
one achieves a well-resolved reconstruction of the 
energy deposited by an interaction. 

Similarly one can remove the charge magnitude 
dependence to reconstruct the interaction depth fl: 

Aq C A + Aq NC A ... 
Aq C A ~ Aq NC A 

The picture of CPG operational principles pre- 
sented in this section can be referred to as zeroth- 
order behavior. First-order effects, which complicate 
this simple picture, are the subject of the following 
sections. 

3. Electron Trapping 

Electron trapping in CZT detectors is a strong 
enough effect that, for COBRA detectors, the signal 
amplitude of an interaction occurring at large depth 
is reduced by about 10% relative to that of an inter- 
action occurring at small depth. This depth depen- 
dence would significantly degrade energy resolution 
were it not corrected. Such a correction is possible 
due to the depth information intrinsic in the raw sig- 
nals. 

P. Luke and E. Eissler introduced a simple weight- 
ing factor (referred to here as w) to modify the dif- 
ference signal in order to correct for electron trap- 
ping Jj: 

Aq C A - wAq NC A (5) 
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Figure 3: CA versus NCA signal heights for 228 Th 
calibration data for one COBRA detector. 7-lines 
from the source appear as upward-sloping lines in 
this plot. The red line indicates 2.615 MeV events 
from 208 T1 decay. The slope of this line, 0.80, is the 
empirically-determined weighting factor w. In the 
absence of trapping this slope would be 1 . 



The weighting factor is empirically determined and 
is always less than one. The correction for trapping 
can be understood qualitatively: as interaction depth 
increases, — Aq^cA becomes smaller while the sig- 
nal reduction due to trapping becomes larger. Under- 
weighting the NCA amplitude thus artificially low- 
ers the amplitudes of low-depth interactions to bet- 
ter match the trapping-degraded amplitudes of high- 
depth interactions. 

The weighting factor w can be determined by an 
optimization procedure seeking the best energy res- 
olution. Alternatively if information from both sig- 
nals is recorded separately the weighting factor can 
be determined from a single calibration run by anal- 
ysis of the CA and NCA signal amplitudes for a cal- 
ibration line. Figure [3] illustrates this procedure as 
performed in COBRA. 

In the following paragraphs an analytical treat- 
ment of the electron trapping effect in CPG detec- 
tors is presented, using a few simplifying idealiza- 
tions. This treatment is similar to the derivation in 



the Appendix of Hfl, although here the anode signals 
are separately derived. First assume a mean trapping 
length A that is valid for the entire electron drift path. 
This is only true to the extent that both the crystal 
properties and the electric field are uniform through- 
out the detector volume. In actuality the electric field 
is nearly uniform for a large majority of the electron 
path for most interaction locations, and thus for most 
of the drift distance in which trapping occurs^ Thus 
the magnitude of the charge as a function of drift 
distance d is taken to be 

Q = Qoe- d/X (6) 

As with zo, A and d are expressed in terms of detec- 
tor length and are thus dimensionless. A is consid- 
ered a free parameter in this treatment. Physically it 
depends upon the electron mobility-lifetime product 
for the crystal and the voltage applied to the crystal, 
both of which vary by detector. Thus A is detector 
dependent. 

A second simplification is to assume that the 
weighting potentials of both anodes take their bulk 
form along the entire electron drift path up until the 
anode plane is reached, at which point the CA and 
NCA weighting potentials abruptly change to 1 and 
0, respectively. The idealized weighting potentials 
as functions of depth z are illustrated in Figure |4] 
In this model the electron drift is reduced to one di- 
mension, with the drift entirely in the negative z di- 
rection. For most interaction locations this is a rea- 
sonable approximation in that, again, most of the 
trapping occurs before deviations in the weighting 
potential become significant and before the electron 
path is diverted by the near-grid field. The expected 
size of the error introduced by this approximation is 
discussed in a later section. 

Given these assumptions the expected anode sig- 
nal amplitudes can be calculated. We divide the cal- 
culation into two parts. The first part corresponds 
to the electrons drifting through the region of uni- 
form weighting potentials, where ^j^- = — 5. Since 
the charge magnitude changes along the drift path 
we calculate this with an integral form of the Ramo 
equation, yielding 



Agi - / -Q e- iz °- z)/x dz 
Jo * 

= iQ A(l-e- z °. 



(7) 



2 This is true considering the electrode geometry if one as- 
sumes that space charge has a negligible effect. The shape of 
charge pulses in COBRA detectors show no evidence of signifi- 
cant field distortions due to space charge. 
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Figure 4: Idealized weighting potentials of CA 
and NCA corresponding to path of charges drifting 
through a CPG detector. 



This result is valid for both the CA and NCA sig- 
nals. It is closely related to the Hecht equation |5, 
p. 489], differing only by the factor one-half. In a 
signal formed by the sum of the CA and NCA sig- 
nals (sometimes referred to as the cathode signal) the 
exact Hecht equation is recovered. This corresponds 
to the weighting potential of a simple parallel-plate 
electrode geometry. 

The second part of the calculation of the anode 
signal amplitudes corresponds to the abrupt changes 
in weighting potential at z=0: +1/2 for CA and -1/2 
for NCA. Inserting the final charge into the Ramo 
equation one finds 



-z /\ 



The full expressions for the anode signals are 



AqcA = \Qo 



AqNCA = 2*5° 



A 1 



-z /X 



-z a /\ 



-z /\ 



(8) 
(9) 

(10) 
(11) 



By inspection we see that, in the limit of large 
A, the corresponding no-trapping forms (Equations [1] 
and|2|i are restored. 

As in the zeroth-order case the position depen- 
dence can be removed by a linear combination of 



the two signals to reconstruct the charge excited into 
the conduction band by the interaction: 



AqcA 



A-l 
A + l 



Aq 



A 



NCA = 



A+l 



Qo 



(12) 



It is noteworthy that the trapping effects, themselves 
of exponential form, are canceled exactly with a lin- 
ear combination. For this it is necessary only that 
the zo-dependent terms in equations [10] and QT| are 
identical in form. On the left-hand side of equa- 
tion [12] is a weighted difference signal equivalent to 
the one introduced by Luke. On the right-hand side 
is the excited charge reduced by a constant factor, as 
expected from the qualitative understanding of the 
weighted difference. The constant is of course ul- 
timately absorbed into the energy calibration. The 
relationship between the parameter A and the empir- 
ical weighting factor w is 



A 



1 



(13) 



1 — w 

An equivalent relationship is derived in Jil. 

It is also possible to construct an analytical ex- 
pression for the interaction depth z$ from the raw 
signal formulas by eliminating Qo from equations 
[10]and[n] 



A In 



1 1 Aq C A + Aq NC A 
A Aq C A - Aq NC A 



(14) 



Note that as expected the zeroth-order depth for- 
mula is recovered in the limit of large A. This 
trapping-corrected depth formula is easily imple- 
mented because the parameter A can be calculated 
from w which is already determined during the en- 
ergy correction procedure. 

4. Evaluation of the trapping model 

In the framework of the simple model for trap- 
ping presented here we can predict the performance 
of the zeroth-order depth formula, now seen as an 
approximation. It is clear that it is non-linear; more- 
over, the zeroth-order formula will overestimate the 
interaction depth in a detector-dependent way. From 
the empirical weighting factors w determined for 
COBRA detectors the overestimation can be ex- 
pected to reach a maximum of around 3-10% for 
zo = 1. 

Figure [5] illustrates the advantages of the new 
depth formula. Sources of alpha radiation at the 
cathode surfaces for the detectors currently in op- 
eration at LNGS produce an excess of events near 
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Figure 5: Interaction depth z of near-cathode events 
in the COBRA data collected at LNGS in 31 detec- 
tors with deposited energies greater than 2 MeV. An 
excess of interactions occurring very near the cath- 
ode is known to exist. Uncorrected depth reconstruc- 
tion results in depths overestimated by a few percent 
and a wider distribution due to detector differences. 



z = 1. The true interaction depth for these events is 
expected to be very close to 1 . The figure compares 
the depth distribution using both the zeroth-order 
and trapping-corrected formulas. The predicted ad- 
vantages are evident: the corrected depths are much 
closer to one, and they are much more tightly dis- 
tributed. 

The empirically determined parameter A can be 
used to estimate the electron mobility-lifetime prod- 
uct (jir) for each detector as follows: 



\xr = 



L 2 X 



(15) 



where L is the detector length (1 cm) and Vbuik 
is calculated to be the potential difference between 
the cathode and the average of the anode potentials 
(HV - ^GB). Figure [6] shows the resulting pa esti- 
mates for 30 COBRA detectors. The mean value of 
1.1 x 10~ 2 cm 2 / V is close to those quoted for CZT 
in several publications (see Table 7 in |9|]). However 
the wide distribution of values suggests significant 
variation across detectors, and possibly additional 
effects unaccounted for by this simplified modelH 




0.03 

lit (cm 2 /V) 

Figure 6: Electron mobility-lifetime products of 
30 COBRA detectors, estimated by the method de- 
scribed in the text. The mean value is 1.1 x 10~ 2 
cm 2 / V. 



5. Additional first-order effects 

5.1. Hole drift 

In CZT holes have much poorer transport prop- 
erties than electrons; nevertheless they contribute 
significantly to the raw signal amplitudes. Because 
holes drift much more slowly than electrons, it is first 
necessary to consider whether their contributions to 
signal amplitudes are fully recorded in the data pro- 
cessing procedure. In standard COBRA pulse pro- 
cessing the effective time window within which the 
amplitude is measured is greater than 3 /zs. The 
lifetime of holes in CZT is estimated at 1 fis, so 
that more than 95% of the hole charge is trapped 
within the processing time window — assuming that 
the holes do not reach an electrode within this time. 
Thus to a reasonable approximation we can assume 
holes contribute fully to the processed signal ampli- 
tudes. 

In analogy to the parameter A in the treatment of 
electron trapping, we take p to be the mean trapping 
length for holes, normalized to the detector length. 
Because typical A values are around 10, and because 
the mobility-lifetime product for holes in CZT is es- 
timated to be about 1% that of electrons, p is esti- 
mated to be around 0. 1 . Thus for many interaction 
depths we do not expect holes to drift far enough to 



3 One potential effect is an imbalance in the CA and NCA 
signal amplification factors arising from the specific properties of 
the electronics. If such an imbalance exists it will contribute to the 
weighting factor w independently of electron transport properties. 



However for this study such imbalances were measured and the A 
and /it values corrected accordingly. Neither can temperature be 
a factor since all detectors are operated under essentially identical 
conditions at room temperature. 
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be collected at an electrode; rather they will be al- 
most completely trapped after drifting a fraction of 
the detector length. 

Qualitatively the hole effect can be pictured as 
follows: holes will drift a short distance toward 
the cathode, along a path in the bulk region of the 
weighting potential. They thus contribute equal pos- 
itive contributions to the CA and NCA signal am- 
plitudes. Since the energy reconstruction is based on 
the difference of these signals, the hole contributions 
will cancel each other. On the other hand the depth 
reconstruction is based on the sum of the anode sig- 
nals, and so the hole effect can lead to a significant 
overestimation of interaction depth. 

Quantitatively we can consider the hole effect us- 
ing the same framework as was used for electron 
trapping. Assuming that the interaction depth occurs 
in the bulk region of the weighting potential, we find 

Wholes = \QoP (l - e-^-^l?) (16) 

This contribution applies to both anode signals. 
Considering the small value of p we see by in- 
spection that the hole contribution is approximately 
\QoP for interaction depths far from the cathode 
(where 1 — zq is several times p). In contrast the 
expression falls to zero as zq approaches 1, since for 
near-cathode interactions the holes have no distance 
to drift. Therefore the hole contribution is zero for 

Zq = 1. 

The contribution of holes to the signals has 
a second-order effect on the energy reconstruc- 
tion, producing a small reduction for near-cathode 
events For interactions far from the cathode 

the weighted difference signal (Equation [TZt is en- 
hanced by approximately p/\, on the order of 1%. 
This constant factor will be absorbed into the energy 
calibration. However the hole enhancement disap- 
pears for interactions very close to the cathode. Thus 
the energy depositon of near-cathode events will be 
underestimated by about 1%. This is comparable to 
the energy resolution of COBRA detectors. The ef- 
fect is noticeable only in a layer whose thickness is 
a few percent of the detector length. Figure [7] shows 
this small effect on a calibration line with roughly 
uniform depth distribution. 

The hole effect on depth reconstruction is larger. 
It can be adequately described by considering the ad- 
dition of the hole contribution to the zeroth-order 
equations (Equation [T] and [U and the correspond- 
ing effect on the zeroth-order depth formula (Equa- 
tion 0J. For interactions far from the cathode, the 
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Figure 7: Distribution of events in depth and energy 
in a COBRA detector irradiated with a 137 Cs source. 
The reconstructed energy of 662 keV 7 events is 
somewhat underestimated at depths near 1 due to the 
hole effect leading to a deformation of the line near 
the top of the plot. 

depth will be overestimated by approximately p. For 
near-cathode interactions there is no overestimation, 
as is clear from Figure [5] Figure [8] illustrates the 
relationship between true and measured depth. The 
clearest evidence of this effect from data is the small 
number of events below a certain measured depth for 
populations of interactions known to have uniform 
depth distributions, as exemplified in Figure [9] 

In principle of course one can also correct for the 
hole effect in the depth reconstruction. However it 
is more difficult to calculate the parameter p in the 
calibration process. Another difficulty is that the 
analytic expression for hole-corrected depth has no 
solution if the zeroth-order depth is greater than 1 . 
Since the finite resolution of measured depth some- 
times results in values greater than 1, this is problem- 
atic for an event-by-event hole correction algorithm. 

5.2. Near-anode distortions 

CPG event reconstruction depends on the nearly 
equal CA and NCA weighting potentials at the inter- 
action point. However the two weighting potentials 
differ significantly at locations near the anode plane, 
as illustrated in Figure [2] 

The corresponding effects are somewhat compli- 
cated, and depend not only on depth zq but also 
on the x and y coordinates of the interaction point. 
The energy can be either underestimated or overes- 
timated depending on whether the interaction point 
lies above a NCA or a CA rail. The thickness of the 
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Figure 8: Relationship between measured depth and 
true depth (red) based on an analytical treatment of 
the hole effect, assuming p = 0.1. The measured 
depth is corrected for electron trapping but not cor- 
rected for the hole effect. 



layer near the anode in which these effects are sig- 
nificant is greater near the lateral surface of the de- 
tected than it is near the center of the anode plane, 
where it is smaller than 0.5 mm. Where it is thinnest, 
holes can actually cancel the distorting effect, by 
drifting toward the cathode into the bulk region of 
the weighting potential. 

Where the effective near-anode layer is thickest, 
holes will not drift far enough to cancel the distorting 
effects. As with the reconstructed energy, the recon- 
structed depth is also distorted, but it will generally 
be smaller than p, so that a judicious depth selection 
can remove such events from consideration in data 
analysis. 

Very near the anode the electric field in the de- 
tector is dominated by the grid bias, which pro- 
duces a different effect. The region of effect depends 
strongly on the applied bias, but is estimated to be 
typically around 0. 1 mm. For interactions in this re- 
gion the holes will drift toward the NCA instead of 
the cathode. The corresponding hole contributions 
to the signals will be opposite in sign, thus strongly 




0.2 0.4 0.6 0.8 1 



Figure 9: Depth distribution of events from a 
COBRA detector for data taken at LNGS in an en- 
ergy range of 180-300 keV. The events consist al- 
most entirely of intrinsic /3-decay of 113 Cd, and 
thus are known to be evenly distributed throughout 
the detector (although some events near the detec- 
tor surfaces will be lost due to escaping f3 particles). 
The low count of measured depths below about 0. 1 
is evidence of the hole effect. The hole effect is 
also partially responsible for the nonuniform distri- 
bution at greater depths. Distortions in energy- and 
depth-reconstruction at lower depths due to the de- 
tails of the weighting potential also contribute to this 
nonuniformity. 



4 The distortions in the weighting potential can be reduced if 
a bias is applied to the guard ring [ 10]. However for simplicity 
COBRA detectors are currently operated with the guard ring un- 
connected. 
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enhancing the difference signal. Due to the short 
hole drift distance, most holes are collected before 
they are trapped. The shape of the weighting poten- 
tials corresponding to the paths of both electrons and 
holes differ significantly from those corresponding 
to events in the bulk. The net effect is that the recon- 
structed energy will be approximately doubled fl. 
The reconstructed depth will be very close to zero. 

Signals from interactions far from the anode are 
also affected by the detailed form of the weighting 
potentials near the anode. If such effects are consid- 
ered Equations [10] and QT| must be corrected. This 
correction will vary with the x and y positions of the 
interaction. However we can approximate the size 
of the correction by considering a specific case: an 
interaction occurring on a line leading to the center 
of a CA rail in the central region of the anode plane. 
In this case the weighting potential as a function of 
z fits well to an exponential form with characteristic 
length a of approximately 0.015 (0.15 mm). We find 
then that Equations[lO]and[TT|are corrected by terms 
of order a/ A. This is a third-order effect, much less 
than 1 %, and thus wholly negligible considering the 
energy and depth resolutions achievable with the de- 
tectors. 

6. Summary 

It has been shown that a simple model for elec- 
tron trapping with an idealized form for the weight- 
ing potential provides a mathematical framework to 
explain the use of a weighted difference signal to re- 
construct event energies in a coplanar grid detector. 
Using the same model a new formula for the inter- 
action depth which corrects for electron trapping ef- 
fects has been introduced. This formula demonstra- 
bly improves the depth reconstruction, especially for 
events near the cathode. Other effects that can dis- 
tort the energy and depth calculations have been dis- 
cussed in an effort to understand event reconstruc- 
tion in CPG detectors as accurately as possible. 
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